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Abstract 

In September 1997 the National Transportation 
Safety Board (NTSB) requested assistance from 
the NASA Lewis Research Center (LeRC) Icing 
Branch in the investigation of an aircraft accident 
that was suspected of being caused by ice 
contamination. In response to the request NASA 
agreed to perform an experimental and 
computational study. 

The main activities that NASA performed were 
LeRC Icing Research Tunnel (IRT) testing to 
define ice shapes and 2-D Navier-Stokes analysis 
to determine the performance degradation that 
those ice shapes would have caused. 

An IRT test was conducted in January 1998. 

Most conditions for the test were based upon raw 
and derived data from the Flight Data Recorder 
(FDR) recovered from the accident and upon the 
current understanding of the meteorological 
conditions near the accident. 


Using a two-dimensional Navier-Stokes code, 
the flow field and resultant lift and drag were 
calculated for the wing section with various ice 
shapes accreted in the IRT test. Before the final 
calculations could be performed extensive 
examinations of geometry smoothing and 
turbulence modeling were conducted. 

The most significant finding of this effort is that 
several of the five-minute ice accretions 
generated in the IRT were found by the Navier- 
Stokes analysis to produce severe lift and drag 
degradation. The information generated by this 
study suggests a possible scenario for the kind of 
control upset recorded in the accident. 

Secondary findings were that the ice shapes 
accreted in the IRT were mostly limited to the 
protected pneumatic boot region of the wing and 
that during testing, activation of the pneumatic 
boots cleared most of the ice. 
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Background 

In early 1 997, a twin-engine commuter aircraft 
(Figure 1 ) crashed during a rapid descent after an 
uncommanded roll excursion. There were no 
survivors. The ground impact forces and post 
accident fire destroyed the aircraft. Instrument 
meteorological conditions prevailed at the time 
of the accident 1 . During the investigation, the 
National Transportation Safety Board (NTSB) 
contacted the Icing Branch at NASA Lewis 
Research Center. The NTSB requested NASA's 
participation since ice contamination was 
suspected as being a factor. NASA offered to 
pursue a study that included experimental and 
computational elements to identify the broad 
nature of possible ice contamination and the 
likely associated perf ormance degradation. 

Project Methodology 

The primary goal of the NASA activity 
supporting the NTSB investigation was to 
attempt to identify a mechanism explaining the 
stability and control upset that preceded the 
accident. The methodology that was adopted to 
satisfy this goal was a multistep process that 
included both computational and experimental 
tasks. 

The first step was to perform initial ice accretion 
predictions with the NASA LEWICE program to 
define the impingement limits and total mass of 
ice to be expected. Since the specific 
meteorological conditions at the time of the 
accident are unknown, the LEWICE calculations 
allowed the generation of ice shapes over the 
range of possible conditions. 

The LEWICE calculations were followed by a 
test in the NASA Lewis Research Center's Icing 
Research Tunnel (IRT). The IRT testing was 
considered necessary since LEWICE, as with all 
ice accretion codes, does not model all of the 
physics of the accretion process. The complex 
movement of liquid water on the airfoil prior to 
its freezing is currently modeled by simplified 
relationships in all ice accretion codes. Due to 
this limitation and other assumptions required by 
ice accretion codes, the ice shapes predicted at 
temperatures near freezing are not always as 
accurate as those predicted for lower 
temperatures. 


The contaminated airfoil geometries measured 
during the IRT test were then analyzed with a 
two-dimensional (2-D) Navier-Stokes code to 
predict their degraded aeroperformance. Due to 
the time constraints of this project, the icing 
wind tunnel test concentrated on measuring the 
ice accretion shapes, with all aerodynamic 
analysis performed computationally. 

To verify that the 2-D analysis was producing 
realist e results (that no significant three- 
dimensional (3-D) flow effects were present), a 
limited 3-D Navier-Stokes analysis was 
conducted. The 3-D analysis was focused on 
assuring that no significant span-wise How was 
present that would invalidate the chord-wise 
flow assumption required by the 2-D Navier- 
Stokes analysis. 


IRT Test 

The Icing Research Tunnel (Figure 2) is a 
closed-loop atmospheric tunnel with a test 
section that is 6 ft high, 9 ft wide, and 20 ft long. 
It is equipped to support testing at airspeeds from 
50 to 430 mph in a temperature and water- 
droplec environment that simulates natural icing 
conditions 2 . Temperature is controllable 
between -20 and +33°F. Clouds can be 
produced in a controlled manner with droplet 
mean volumetric diameter (MVD) between 9.5 
and 270 pm. 

For this test an existing model was utilized. The 
model was a modified piece of wing from the 
subjec aircraft centered on the inboard region of 
the aileron. The model had a 6-foot span and 
was mounted vertically in the IRT (Figure 3). 

The wmg had a modified trailing edge and 
pneumatic boot. The trailing edge was faired 
behind the aft spar to avoid the mechanical 
complexity of the aileron. The pneumatic boot 
was modified for its original test to allow the 
incorp oration of thermal measurements. This 
featurt was not utilized for this subject testing, 
however, the modification resulted in the boot 
extending further aft on the lower surface than a 
produc tion boot for this wing. This extended 
boot surface acted as an enhanced ice accretion 
location that is evident when reviewing the test 
ice sh^pe tracings. Except for the cases when the 
pneumatic boot was activated, the manufacturers 
recommended level of vacuum was applied to 
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the boot throughout the testing to prevent auto- 
inflation. To ensure that the air temperature in 
the test section was accurate at the desired near- 
freezing condition, an auxiliary temperature 
probe was mounted near the model and is visible 
in Figure 3. 

Test Environment Conditions 

Based upon the findings of the investigation’s 
meteorological team*, the accretion time, the 
maximum liquid water content (LWC), and the 
cloud droplet’s median volumetric diameter 
(MVD) were estimated for this test. The 
accretion time was estimated at 5 minutes. This 
was the length of time required by the accident 
aircraft to transit an area of low reflectivity 
evident on regional weather radar. This time 
also agrees with estimates based upon a 
measurable drag increase calculated by the 
investigation’s performance team from flight 
data recorder (FDR) information 3 4 . The 
maximum LWC was estimated to be 0.8 g/nf 
based upon microwave radiometer measurements 
of the same airmass following the accident. The 
maximum MVD of the cloud droplets was 
estimated to be no greater than 70 to 80 pm. 

This droplet-sizing estimate was considered 
reasonable since the aircraft was known to have 
been in the clouds, and droplets larger than that 
are normally associated with below-cloud 
drizzle 5 . 

The remaining test parameters were defined 
based upon a combination of the accident 
meteorological data and the FDR inlormation. 
The times that defined the assumed icing period 
from the meteorological report were used to 
access the pertinent information from the FDR. 
Representative airspeed, aircraft angle-of-attack 
and air temperature values were calculated by 
averaging over the assumed ice accretion period 
(total air temperature was derived from static air 
temperature measurements and corrected tor 
altitude and airspeed, and wing angle-of-attack 
was calculated by adding the wing incident angle 
to the aircraft body angle-of-attack). 

The combination of the assumed cloud 
conditions and the measured flight environment 
data was used to define the IRT test matrix 
(Table 1). Four conditions were considered the 
baseline series in the matrix. Due to 


uncertainties in the FDR measured values of 
angle-of-attack and temperature, additional series 
were added to the matrix that examined higher 
and lower angle-of-attack than the baseline value 
and higher and lower temperatures than that 
baseline value. An additional series with longer 
accretion time was also added since the ice 
accretion time was inferred information and not 
directly measured. A series of supercooled large 
droplet (SLD) conditions was added during the 
testing at the request of the NTSB to ensure the 
completeness of the database. To add to the 
level of confidence in the IRT results, all primary 
conditions were repeated at least once. 

Test Procedure 

The IRT test was conducted over two weeks in 
January 1998. The conditions described in Table 
1 were used for this test and the most significant 
conditions were repeated at least once, as seen in 
Table 2. The normal process used during this 
test was to cool the IRT to the desired 
temperature and set the airspeed. When the 
temperature and airspeed were stable at the 
desired test conditions, the spray was activated 
for the length of time at air and water pressures 
defined by the test matrix to simulate a natural 
icing cloud. When the spray was shut off, 
several wake survey probe traverses were made 
behind the model. These wake probe traverses 
were used to calculate an uncorrected tunnel 
centerline section drag. Due to project time 
constraints, further correction, analysis, and 
comparisons of this data have not been 
completed. Following the wake probe traverses, 
the tunnel fan was stopped so that test personnel 
could enter the test section. Upon entering the 
test section, the test engineers would first 
photograph the model. Both digital and film 
cameras were utilized for this test, with the film 
photos typically being used for close-up detailed 
shots and the digital photos for overview shots 
(Figures 4 and 5). When the photographs were 
completed, the leading edge ice accretion was 
sliced with a hot knife, then traced onto a 
cardboard template (Figure 6). These tracings 
were made at 30, 36, and 42 inches from the 
wind tunnel floor (the 36 inch location is the 
vertical center of the test section). The ice 
shapes originally traced onto cardboard are then 
later digitized and stored electronically. All 
tracings presented here are digitized data 
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(Figure 7). Three channels of video were also 
recorded for each ice accretion. The video 
cameras captured upper and lower leading edge 
views along with a broad model view. The video 
signals were recorded to SVHS format tapes. 

Resultant Ice Accretions 

Repeatability was good for all cases, with the 
volume of ice, impingement limits, and horn or 
ridge location closely matched for repeat cases. 
The majority of the ice shapes accreted during 
this test fell into two categories. The first was a 
thin ice shape with a rough surface made up of 
small bumps of varying size and sharpness 
(Figure 8). The second category of ice shape had 
the same general volume of ice with similar 
roughness elements as the first category, but also 
a prominent ridge (Figure 9). This ridge was 
approximately 0.5 inches high. 

The following is a description of the trends 
apparent in the test data: 

• The baseline series (see Table I) resulted in 
thin ice accretions with distributed 
roughness and no ridge formation for all 
four of the series conditions (Figures 10 to 
13). 

• Compared to the baseline series, the higher 
angle-of-attack series (at 7 degrees instead 
of 5 degrees) produced ice shapes that were 
shifted towards the lower surface, as would 
be expected with this change of angle-of- 
attack (Figure 14). Since the upper 
impingement limit moved forward, the 
overall ice thickness of these shapes was 
somewhat greater than that seen in the 
baseline series. 

• The lower angle-of-attack series cases (at 3 
degrees instead of 5 degrees) were not 
significantly different than the baseline 
series. 

• The higher temperature series cases (3 1 
degrees Fahrenheit total temperature instead 
of 30 degrees) were not significantly 
different than the baseline series, but may 
have had a very slight decrease in the total 
volume of ice accreted. 


The lower temperature series (at 28 degrees 
Fahrenheit total temperature instead of 30 
degrees) produced slightly more total ice 
volume than the baselines series and more 
importantly started to exhibit some ridge 
formation for the 40 and 70 pm MVD cases 
(Figure 15 and 16). 

The second lower temperature series (at 26 
degrees Fahrenheit total temperature instead 
of 30 degrees) exhibited an ice ridge for all 
four conditions (Figures 17 to 20). The 
ridges apparent in these accretions are 
different from a standard horn formation in 
that they are surrounded by only a thin layer 
ol ice, where a standard horn is part of a 
more substantial ice accretion. The ridge 
w is most prominent for the two cases with 
LWC=0.8 g/m 3 (matrix numbers 26 and 27, 
Figure 17 and 18). For the two cases with 
constant LWC, the ridge moved aft with 
increasing MVD. For the two cases with 
constant MVD, the ridge moved aft with 
increasing LWC. 

In the longer time series (with an accretion 
time of 10 minutes instead of 5), self- 
shedding occurred in matrix number 32 
(Figure 21). For the remaining cases in this 
series the ice mass was greater than in the 
baseline series, as would be expected. 

These shapes also exhibited a ridge 
formation (Figures 22 to 24). However 
th ;se ridge formations appeared to be much 
more like a standard horn formation. These 
ice shapes were in general larger and 
rougher than those seen in the baseline 
series. 

T1 e additional supercooled large droplet 
se ies (which extended the MVD up to 270 
pr i from the baseline maximum of 70) 
sh iwed no significant difference in the 
leading edge ice accretion over those 
exhibited in the baseline series, except that 
the 175 and 270 pm cases self-shed. 

Ti e pneumatic boot was operated for several 
ca ;es in the de-ice series. It was operated for 
the matrix number 3, 26, and 27 conditions. 
Depending on the case, the pneumatic boot 
was either activated at spray-off only 
(Figure 25), at two minutes into a five- 
minute spray (Figures 26 and 27), or at 
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two minutes and at spray -off (Figures 28 
and 29). In all cases when the pneumatic 
boot was operated at the end of the icing 
spray, the leading edge ice shed very well. 
For the two cases when the boot was 
activated at two minutes into a five-minute 
spray, in one case it prevented the 
development of a ridge formation, and in the 
other case it did not. Whether the ridge was 
suppressed or not was probably dependant 
on how well the boot cleared the ice at the 
two-minute activation. It should be noted 
that the pneumatic boot used during this test 
was not a production boot and it had never 
experienced the operational environment. 
However, no special steps were taken to 
enhance the boot performance, such as the 
application of Ice-X, and it was operated at 
the same pressures as a flight system. The 
boot performance in the deicing tests is 
therefore considered to be broadly 
representative of that of a flight system. 

The physics involved in the development of the 
ice ridges observed during this test are not well 
understood. It is very possible that the ridges are 
a result of the damming of surface runback water 
by ice feathers. The ridges appear to have 
formed at the aft edge of the primary ice 
accretion where ice feathers begin. Ice damming 
by ice feathers has been observed by close-up 
videography 5 6 , and may explain this form of ice 
growth. More research is required into the 
formation of the ice ridges seen in this test and, 
more generally, into the movement of liquid 
water on developing ice accretions. 

Computational Analysis 

The computational analysis performed as part of 
this effort is described in two A1AA reports™. 
The computational analysis consisted of both 2- 
D and 3-D Navier-Stokes 9 computations. As 
described above, the performance of the ice- 
contaminated wing was calculated via 
computational techniques instead of by 
experiment in order to complete the analysis in a 
timely manner for the NTSB investigation. 

Four cases of the IRT ice shapes were selected 
for computational analysis based upon their 
varying nature. The ice shape from test matrix 
condition number 2 was selected since it 


exhibited a thin ice accretion with many small 
roughness elements. Matrix number 3 was also 
selected since it exhibited larger, sharper 
roughness elements than did matrix number 2. 
Matrix number 9 was selected because it 
exhibited the roughness of matrix number 2, but 
also had a larger accretion on the upper surface. 
Finally, the resultant ice shape form matrix 
number 26 was selected since it exhibited a 
significant ridge formation. 

For the 2-D calculations, a significant effort was 
undertaken to determine the best level of ice 
tracing point smoothing and optimal grid point 
resolution™. The calculated lift curves of the 
four ice shapes examined compared to that of the 
clean airfoil are presented in Figure 30. As 
expected, when ice was added to the airfoil, the 
calculated lift was degraded. The matrix 2 shape 
had the least degradation since it had the lowest 
level of roughness of the shapes examined, 
followed by the matrix 3 shape, which had 
larger, rougher roughness elements. When the 
ice thickness increased, as in matrix number 9, 
the lift degradation was further increased. The 
surprising result was the level of degradation 
when a ridged ice shape was examined. The 
major difference between the ice shape from 
matrix number 2 and that from matrix number 26 
was the ice ridge. The roughness levels observed 
in the two cases were very close. When the ice 
ridge was present, a 45 percent reduction of 
maximum lift coefficient (C, max ) was predicted 
in comparison to the clean airfoil. Since the 
ridged ice shape caused the greatest lift 
degradation, it was selected for all further 
computational analysis. 

Of equal significance compared to the lift 
degradation caused by the ice contamination is 
the effect of aileron deflection for contaminated 
airfoils. Navier-Stokes calculations were 
performed for the wing in the configuration at 
the moment of the control upset. The moment of 
control upset was defined for this activity as the 
point when the aircraft’s roll no longer tracked 
the aileron input. The roll angle and aileron 
(control wheel) position were recorded on the 
FDR (Figure 31.). The configuration of the 
aircraft with regards to bank and aileron position 
is shown in Figure 32. It should be noted that 
the aircraft flight variables of interest here 
(airspeed, angle-of-attack, altitude, and 
temperature) are different from those used to 
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simulate the possible ice shapes in the IRT, since 
the accretion of ice and the control upset 
occurred at different times and spatial position. 

When the aileron detlection effect was added to 
the degradation caused by the ridged ice shape 
the lift curves in Figure 33 resulted. The clean 
left wing (Clean LW) line represents the lift 
curve for the airfoil with the aileron deflected 
down at 2.56 degrees. The clean right wing 
(Clean RW) line represents the lift curve for the 
airfoil with the aileron deflected up at 2.74 
degrees. The left wing generates more lift than 
the right; therefore the aircraft will roll right. 

This was the intention of the pilot when these 
control inputs were made. However, with the 
ridged ice shape contamination on the wings, 
both the Iced LW and Iced RW lines show how 
the left and right airfoil lift curves were 
degraded. As angle-of-attack increases the 
differential lift decreases and at an angle-of- 
attack around 1 1 degrees the difference in lift 
between the two wings goes to zero. This 
degraded roll control authority is a possible 
answer as to why the aircraft continued to roll to 
the left when the aileron inputs were 
commanding a right roll. 

At the completion of the 2-D analysis a limited 
3-D analysis was undertaken 8 . The left wing 
geometry, including aileron deflection but minus 
engine nacelle, was applied to a 3-D structured 
grid. Developing a grid that accurately captured 
the wing geometry and also allowed the Navier- 
vStokes code to converge took a great deal of 
effort. To allow the convergence of the Navier- 
Stokes code in a timely manner the grid 
resolution used was relaxed. The resultant 3-D 
flow field demonstrates a chord- wise flow over 
the wing. The trailing edge stall mechanism is 
also similar to that seen in the 2-D analysis, 
however, the stall occurred much more abruptly 
in 3-D than it did in 2-D. In 3-D the separation 
started right before the stall, but for 2-D the 
separation started several degrees before the 
stall. The angle-of-attack for stall in the 3-D 
analysis is also higher than that for the 2-D 
analysis. These differences may be caused by the 
relaxed grid resolution utilized for the 3-D 
calculations. While explainable, the difference 
between the 2-D and 3-D results does warrant 
further investigation. A preliminary 
experimental comparison has also been made 
with results generated by the University Of 


Illinoi: at Urbana/Champaign 10 . While the 
airfoil used for the experiment was different than 
that used for this computational activity, the ice 
ridged shape was similar. The experimental stall 
angle is greater than that predicted by the NASA 
2-D ef fort, but occurs due to trailing edge flow 7 
separation. The nature of the flow results appear 
to be similar to the NASA results, but as stated 
before more research will be required before the 
computational results can be utilized 
independently with complete confidence. 

Summary of Significant Experimental and 
Computational Finding s 

The ice shapes accreted during the IRT test 
largely fell into two categories. The first was a 
thin ies shape with a rough surface made up of 
small oumps of varying size and sharpness. The 
second category of ice shape had the same 
general volume of ice with similar roughness 
elements as the first category, but also a 
prominent ridge. This ridge was approximately 
0.5 inc hes high. 

In all cases when the pneumatic bool was 
operated at the end of the icing spray, the leading 
edge i ;e shed very well. 

Of the four ice shapes analyzed with the Navier- 
Stokes code, the degraded C, max from the ice 
shape with the lowest level of roughness, matrix 
number 2, was about 75 percent of the clean 
airfoil The C, max from the ice shape with larger, 
sharpt r roughness elements, matrix 3, was about 
68 peicent of the clean airfoil. The C, max from 
the ro igh shape with more ice accretion on the 
upper surface of matrix number 9 was about 60 
percent of the clean airfoil. And the C t max from 
the ice ridge case, matrix 26, was about 53 
percent of the clean airfoil. 

With *he aileron effects added to the lift 
degrai lation caused by the ridged ice shape, the 
result; tnt lift curves for the left and right wing 
airfoils display that little or no lift difference 
may h ave been present and therefore roll control 
autho ity may have been significantly degraded. 

While 3-D computational and experimental 
result . do not directly confirm the 2-D Navier- 
Stoke ; results, they do confirm the general flow 
characteristics, and stall mechanism. And 
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enough differences exist between the 2-D and 3- 
D grids and the computational and experimental 
geometries to account for the differences in 
results. 


Concluding Remarks 

There are many lessons to be learned by this 
activity. The NASA Icing Branch has gained a 
significant amount of experience in measuring 
wind tunnel generated ice shapes and 
determining the aerodynamic performance for a 
contaminated wing. Along with this added 
knowledge have come more questions. What are 
the physical processes that led to the generation 
of the ridged ice shapes? How do varying ridge 
locations and heights influence airfoil and wing 
performance? How well do the Navier-Stokes 
codes compare to experimental aeroperformance 
data, and how can the agreement be improved? 
And, what is the best way of using the 2-D and 
3-D codes together? There is a significant 
amount of work that needs to be accomplished 
before these question are answered and before 
we have a good understanding and confidence in 
all the methodologies. 

The authors therefore recommend that the 
following activities be pursued: 

• Additional testing be conducted to gain an 
understanding of the ice ridge formation 
process. 

• More analysis and testing be conducted to 
understand the criticality of the ice ridge in 
regards to airfoil performance. This needs 
to include both ridge location and height. 

• Compare the computational results from 2-D 
and 3-D Navier-Stokes code utilized to other 
Navier-Stokes codes using the same 
geometry and conditions. 

• Aeroperformance testing and comparison to 
the Navier-Stokes solutions. 

The other lessons learned through this activity 
relate to the operation of aircraft in icing 
conditions. Based upon the cockpit voice 
recorder from the accident aircraft, it appears 
that the flight crew was unaware of the potential 
hazard of their flight situation. While the total 
ice accretion was probably quite small, the 
aerodynamic effects appear to have been severe. 
For this incident, relying on past experiences and 


guidelines that suggested waiting until 
significant accretions existed before activating 
the pneumatic ice protection systems may have 
been a fatal error. Due to the limited information 
recorded, we will never be completely sure of the 
cause of this accident and can surely never know 
the level of awareness and intentions of the flight 
crew. However, based upon a review of the 
NASA results from this study along with the 
group consensus at a recent workshop 11 , two 
operational suggestions become apparent: 

• Even ice accretions that may appear small 
and benign can be truly hazardous. Treat all 
icing conditions and resultant ice accretions 
as threats to your aircraft. 

• Unless otherwise instructed by the aircraft 
flight manual, activate the ice protection 
systems at first indication of ice accretion; 
do not wait until a significant accretion has 
formed. 
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Tables 


Table 1, Planned Test Matrix 


Matrix # 

A/S 

(knots) 

Total 

Temp 

(F) 

AOA 

(deg) 

LWC 

(g/m’) 

MVD 

(|*m) 

Pair 

(psia) 

DP 

(psid) 

Time 

(min) 

Comments 

1 

172 








Clean Wake Run 

2 

172 

30 

5 

0.8 

20 

37 

98 

5 

Baseline series 

3 

172 

30 

5 

0.8 

40 

19 

83 

5 

Baseline series 

4 

172 

30 

5 

0.52 

40 

10 


5 

Baseline series 

5 

172 

30 

5 

0.58 

70 

8.2 

32 

5 

Baseline series 

6 

172 

30 

7 

0.8 

20 

37 

98 

5 

Higher AOA series 

7 

172 

30 

7 

0.8 

40 

19 


5 

Higher AOA series 

8 

172 

30 

7 

0.52 

40 


33 

5 

Higher AOA series 

9 

172 

30 

7 

0.58 

70 

8.2 

32 

5 

Higher AOA series 

10 

172 








Clean Wake Run 

1 1 

172 








Clean Wake Run 

12 

172 

30 

3 

0.8 

20 

37 

98 

5 

Lower AOA series 

13 

172 

30 

3 

0.8 

40 



5 

Lower AOA series 

14 

172 

30 

3 

0.52 

40 


33 

5 

Lower AOA series 

15 


30 

3 

0.58 

70 


32 

5 

Lower AOA series 

16 

■ 

31 

5 

0.8 

20 


98 

5 

Higher Temp series 

17 

172 

31 

5 

0.8 

40 

19 


5 

Higher Temp series 

18 

172 

31 

5 

0.52 

40 



5 

Higher Temp series 

19 

172 

31 

5 

0.58 

70 

8.2 

32 

5 

Higher Temp series 

■■ 

172 








Clean Wake Run 

■n 

172 








Clean Wake Run 


172 

28 

5 

0.8 

20 

37 

— 

5 

Lower Temp series 

HS 

172 

28 

5 

0.8 

40 

19 


5 

Lower Temp series 

24 

172 

28 

5 

0.52 

40 

30 

33 

5 

Lower Temp series 

25 

172 

28 

5 

0.58 

70 

8.2 

32 

5 

Lower Temp series 

BIB 

172 

26 

5 

0.8 

20 


98 

5 

Lower Temp series II 


172 

26 

5 

0.8 

40 

19 

83 

5 


Krm 

172 

26 

5 

0.52 

40 


33 

5 



172 

26 

5 

0.58 

70 

8.2 

32 

5 

Lower Temp series 11 

■ 

172 








Clean Wake Run 

■n 

172 








Clean Wake Run 

32 

172 

30 

5 

0.8 

20 

37 

98 

10 

Longer Time series 


172 

30 

5 

0.8 

40 


IBS 


Longer Time series 

IBi 

IBB 

30 

5 

0.52 

40 


33 


Longer Time series 

35 

\msam 

30 

5 

0.58 

70 

8.2 

32 

10 

Longer Time series 

— 

172 








Clean Wake Run 


172 

30 

5 

0.6 

100 

6 

BS 

5 

Additional SLD series 

BEK 

172 


5 

0.6 

120 

5 

1 28 

5 

Additional SLD series 


172 

30 

5 

0.85 

175 

5 


5 

Additional SLD series 

■ 

172 

30 

5 

0.85 

270 

2 

1 22 

5 

Additional SLD series 

■n 

172 

28 

5 

0.6 

100 

6 


5 

Additional SLD series 

■m 

172 

30 

3 

0.6 

100 

6 

mm 

5 

Additional SLD series 

■ 

172 

30 

7 

0.6 

100 

6 

| 30 

5 

Additional SLD series 
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Table 2, As-Run Test Matrix 


Matrix # 

Comments 

Date 

Run# 

Repeat 

Date 

Repeat 
Run # 

Repeat#2 

Date 

Repeat#2 
Run # 

1 

Clean Wake Run 

1/20/98 

1 





2 

Baseline series 


2 

1/26/98 

8 

1/26/98 

9 

3 

Baseline series 

1/20/98 

3 


6 

1/26/98 

7 

4 

Baseline series 

1/20/98 

4 

1/26/98 

4 

1/26/98 

5 

5 

Baseline series 

1/20/98 

5 

1/26/98 

2 

1/26/98 

3 

6 

Higher AOA series 

1/21/98 

2 

1/29/98 

7 



7 

Higher AOA series 

1/21/98 

3 

1/29/98 

8 



8 

Higher AOA series 

1/21/98 

4 

1/30/98 

2 



9 

Higher AOA series 

1/21/98 

5 

1/30/98 

3 



12 

Lower AOA series 


6 

1/29/98 

2 



13 

Lower AOA series 

1/21/98 

7 

1/29/98 

3 



14 

Lower AOA series 


8 

1/29/98 

4 



15 

Lower AOA series 


2 

1/22/98 

3 

1/29/98 

5 

16 

Higher Temp series 


4 





17 

Higher Temp series 


5 





18 

Higher Temp series 


6 





19 

Higher Temp series 

1/22/98 

7 





22 

Lower Temp series 

1/22/98 

8 

1/28/98 

2 



23 

Lower Temp series 


9 

1/28/98 

3 



24 

Lower Temp series 

1/23/98 

2 

1/28/98 

4 



25 

Lower Temp series 


3 

1/28/98 

5 



26 

Lower Temp series II 

1/23/98 

4 

1/27/98 

6 




Lower Temp series II 

1/23/98 

5 

1/28/98 

9 




Lower Temp series II 

1/23/98 

6 

1/28/98 

7 



29 

Lower Temp series II 

1/23/98 

7 

1/28/98 

8 



32 

Longer Time series 

1/23/98 

8 





33 

Longer Time series 


9 





34 

Longer Time series 

1/23/98 

10 





35 

Longer Time series 


11 





37 

Additional SLD series 


2 





38 

Additional SLD series 


3 





39 

Additional SLD series 

1/27/98 

4 





40 

Additional SLD series 


5 





41 

Additional SLD series 

1/28/98 

6 





42 

Additional SLD series 

1/29/98 

6 





43 

Additional SLD series 

1/30/98 

4 





3 Deice 

Matrix #3 

With De-ice @ 5min 

1/26/98 

10 





26Deice 

Matrix #26 
With De-ice <s> 2min 

1/27/98 

7 





27 Deice 

Matrix #27 
With De-ice @ 2min 

1/27/98 

8 





27 Deice 

Matrix #27 

With De-ice @ 2min and 5min 

1/27/98 






26Deice 

Matrix # 26 

With De-ice @ 2min and 5min 

1/27/98 

10 






10 

American Institute of Aeronautics an i Astronautics 

































































































11 

American Institute of Aeronautics and Astronautics 





Figure 6., Ice Tracing Technique 



Figure 7., Example Digitized Ice Shape 
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^igure 9., Tracing of Ridged Ice Shape 



Figure 10., Baseline condition, 
Matrix number 2 tracing 



Figure 1 1 Baseline condition, 
Matrix number 3 tracing 


Figure 8., Tracing of Non-Ridge Shape 
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Figure 12., Baseline condition. 
Matrix number 4 tracing 


Figure 15., Lower temperature condition. 
Matrix number 23 tracing 
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Figure 13., Baseline condition. 
Matrix number 5 tracing 


Figure 16., Lower temperature condition, 
Matrix number 25 tracing 



Figure 14., Higher AOA condition. 
Matrix number 6 tracing 


Figure 17., Second lower temperature condition. 
Matrix number 26 tracing 
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Figure 18., Second lower temperature condition. 
Matrix number 27 tracing 


Figure 21., Ten minute accretion condition. 
Matrix number 32 tracing 



Figure 19., Second lower temperature condition, 
Matrix number 28 tracing 


Figure 22., Ten minute accretion condition. 
Matrix number 33 tracing 
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Figure 20., Second lower temperature condition, Figure 23., Ten minute accretion condition. 

Matrix number 29 tracing Matrix number 34 tracing 
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Figure 24., Ten minute accretion condition, 
Matrix number 35 tracing 
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Figure 27., De-ice at 2 minutes of spray. 
Matrix number 27Deice tracing 
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Figure 25., End of spray De-ice, Figure 28., De-ice at 2 minutes and at end of 

Matrix number 3Deice tracing spray, Matrix number 27Deice tracing 
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Figure 26., De-ice at 2 minutes of spray. 
Matrix number 26Deice tracing 
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Figure 29., De-ice at 2 minutes and at end of 
spray. Matrix number 26Deice tracing 
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Figure 30., Lift curves for clean airfoil and 
airfoil contaminated with various ice shapes 


Wh«ol Position- Roll Angle 



Figure 3 1 Flight Data Recorder Tracing 
including upset 



F igure 32., Aircraft and aileron position 
(not to scale) at upset (view from behind 
looking forward) 



AOA idagrses) 


Fij ure 33., Lift curves for clean airfoil and 
ridge ice shape contaminated airfoil 
with aileron deflections 
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